% Version 2.11, April 2022.
%
% Jan Dirk Jansen, Delft University of Technology (TU Delft)
% j.d.jansen@tudelft.nl

These Matlab files have been used to generate the plots for:

Jansen, J.D. and Meulenbroek, B., 2022: "Induced Aseismic Slip and the Onset of Seismicity in Displaced Faults". Submitted for publication to the Netherlands Journal of Geosciences.

The main file is fault_slip_main.m. It contains a set of input parameters to compute and plot depletion-induced stresses and fault slip along an inclined fault for user-defined parameters. 

If parameter flag_paper (in line 107) is set equal to 1, the program produces nearly all the plots from the paper (plus several more). However, the production of Figure 7 is very time consuming and has therefore been switched off. To switch it on again, open file fault_slip_patch_bounds_paper.m and set parameter flag_plot_offset_sensitivity (in line 17) equal to 1. 

Note 1: This is research code which may crash when using different parameters than the ones in the main file.

Note 2: The Matlab files have been developed alongside with the writing of the paper over a period of about two years. During that time the notation in the paper has been sligtly adjusted such that there are small discrepances between the notation in the paper and the code. In particular:
* The fault normal stresses, effective normal stresses and shear stresses are indicated in the paper with the Latex variables \sigma_{\perp}, \sigma'_{\perp} and \sigma_{\parallel} respectively. In the Matlab files they are indicated with sigma_norm, sigma_norm_eff and sigma_shear.
* The initial, incremental and combined stresses are indicated in the paper with the Latex variables \sigma^0_{xxx}, \sigma_{xxx} and \Sigma_{xxx} respectively with xxx representing \perp or \parallel. In the Matlab files they are indicated with sigma_yyy_0, sigma_yyy and sigma_yyy_comb, with yyy representing norm, norm_eff or shear.
* The number of terms in Chebyshev expansions in the paper is taken as N+1. In the Matlab code we use just N terms. For a practically relevant number of terms (typically 100 or more), the difference is neglegible.     
 